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Kopeliovich's suggestion |lj to perform nuclear geometry (Glauber) calculations using different 
cross sections according to the experimental configuration is quite different from the standard prac- 
tice of the last 20 years and leads to a different nuclear geometry definition for each experiment. 
The standard procedure for experimentalists is to perform the nuclear geometry calculation using 
the total inelastic N-N cross section, which results in a common nuclear geometry definition for all 
experiments. The incomplete acceptance of indivdual experiments is taken into account by correct- 
ing the detector response for the probability of measuring zero for an inelastic collision, which can 
often be determined experimentally. 0] This clearly separates experimental issues such as different 
acceptances from theoretical issues which should apply in general to all experiments. Extreme- 
Independent models are used to illustrate the conditions for which the two methods give consistent 
or inconsistent results. 

PACS numbers: 25.75.-q,25.40.-h,24.85.+p,13.85.Hd 



I. INTRODUCTION 

In a recent article, Kopeliovich has discussed 
the case of nuclear geometry calculations for detec- 
tors that are not fully sensitive to the total inelastic 
nucleon-nucleon (N-N) cross section, i.e. those with 
limited acceptance. He suggests that the nuclear 
geometry (Glauber) calculation should be done us- 
ing different inelastic N-N cross sections according 
to whether or not the detector is sensitive to inelas- 
tic N-N difFractive processes. Explicitly, Kopeliovich 
suggests to use the N-N non-diffractive cross section, 
which he estimates to be 30mb, in the nuclear ge- 
ometry calculation for the PHOBOS and PHENIX 
detectors, instead of the 42mb total inelastic N-N 
cross section used by both experiments. 0] The gen- 
eral practice by experimentalists is to use the to- 
tal inelastic N-N cross section in the nuclear ge- 
ometry calculation and then correct for the limited 
phase-space coverage (the acceptance) at the detec- 
tor level 0, IE This practice is based primar- 
ily on the fact that the nuclear geometry calcula- 
tion is then the same for all experiments and so 
can be easily compared. Secondly, in the case of 
fixed target experiments, the number of projectile 
participants, a key nuclear geometry parameter, can 
be directly measured using a Zero Degree Calorime- 
ter, la 13)0 Thirdly, all experiments have acceptance 
effects which must be corrected in any case. 

The issue of non-diffractive versus total inelastic 
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N-N cross section is relevant to the definition of pro- 
jectile participants which can be directly measured 
using Zero Degree Calorimeters (ZDC). The solid 
angle of the ZDC is set to a very small forward cone 
around the beam direction i] ^ y&eam^ 
tect only projectile spectators. The ideal aperture 
would allow the ZDC to measure the full kinetic en- 
ergy of a projectile nucleon in the case of no interac- 
tion, and to measure zero energy for any nucleon in 
the projectile that suffered an inelastic collision, in- 
cluding diffraction dissociation. Any inelastic inter- 
action, including diffraction excitation, which causes 
a projectile nucleon to acquire transverse momen- 
tum and lose energy, moves it out of the ideal ZDC 
aperture. Thus, the energy recorded in the ZDC is 
proportional to the number of non-interacting nu- 
cleons ("spectators") in the projectile, so that the 
number of nucleons which have interacted ("projec- 
tile participants") is straightforwardly deduced. 

It is important to note that at RHIC, the ZDC's 
do not satisfy the ideal criterion and so do not di- 
rectly measure the number of projectile spectators, 
hence participants. Nevertheless, it is preferable to 
keep the nuclear geometry calculations general, in 
case, for instance, improved ZDC's are installed at 
RHIC some time in the future. Furthermore, the use 
of a common nuclear geometry definition by all four 
RHIC experiments 2] has made comparisons among 
them straightforward and transparent. 

An additional argument against adjusting the 
cross section used in the nuclear geometry calcula- 
tion to account for the experimental configuration 
is that this procedure does not work in general and 
may, in certain models, give a different answer from 
the standard method of including the acceptance 
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in the detector response function. This will be il- 
lustrated using examples from extreme-independent 
models which offer the advantage of explicitly sep- 
arating the effect of the nuclear geometry from the 
detector response. 



II. EXTREME-INDEPENDENT MODELS 

A. Standard procedure to correct for 
incomplete acceptance 

The beauty and utility of having a generally ap- 
plicable nuclear geometry calculation, rather than 
one tailored for a particular experiment, is that in 
the extreme-independent-collision models of B-fA 
nuclear scattering, such as the Wounded Nucleon 
Model 13, I3 and Wounded Projectile Nucleon 
Model [3, 13] , the effect of the nuclear geometry of 
the interaction can be calculated independently of 
the dynamics of particle production, which can be 
taken directly from experimental measurements. In 
these models, the nuclear geometry is represented 
as the relative probability w„ per interaction for a 
given number n of total participants (WNM), projec- 
tile participants (WPNM), or other basic elements 
of particle production such as wounded projectile 
quarks (Additive Quark Model- AQM) |l3| or binary 
nucleon- nucleon collisions (NCM), integrated over 
the impact parameter of the B-l-A reaction. Typ- 
ically, Woods-Saxon densities are used for both the 
projectile and target nuclei, and the nucleon-nucleon 
inelastic cross section ct;^^ appropriate to the cm. 
energy of the collision is taken. At AGS energies, 
30 mb was used, corresponding to a nucleon- 
nucleon mean free path of ~2.2 fm at nuclear den- 
sity, while at RHIC, an inelastic N-N cross section 
of 42 mb is appropriate for ^/snn = 200 GeV. ||Tl| 
Once the nuclear geometry is specified in this man- 
ner, experimental measurements can be used to de- 
rive the distribution (in the actual detector) of Et 



or multiplicity (or other additive quantity) for the 
elementary collision process, i.e. a wounded nucleon 
(WNM) or a wounded projectile nucleon (WPNM) , 
a wounded projectile quark (AQM) or an N-N colli- 
sion (NCM), which is then used as the basis of the 
analysis of a nuclear scattering as the result of mul- 
tiple independent elementary collision processes. 

To illustrate the effect of the detector acceptance, 
we use the number of collision model (NCM) as an 
example to calculate an Et distribution. The NCM 
calculation for a B-|-A reaction is given by the sum: 

i-^) ^<TBAj2^nPn{ET) (1) 

where aBA is the inelastic B-l-A cross section, w„ 
is the relative probability for n binary collisions in 
the B+A reaction, from 1 to iV,nax, and PniEx) is 
the calculated Et distribution on the detector for n 
independent N-N collisions. If /i(i?T) is the mea- 
sured Et spectrum on the detector for one detected 
N-N collision, and po is the probability for an N-N 
collision to produce no signal in the detector, then, 
the correctly normalized Et distribution for one N- 
N collision is: 

Pi{Et) = {l-po)fiiET)+PoS{ET) , (2) 

where 5{Et) is the Dirac delta function and 
/ Ji{Et) dET = 1. -P„(£'t) (including the po effect) 
is obtained by convoluting Pi{Et) with itself n — 1 
times 

P^(ET)^j:j^^^,pr\l-PoTMET) , (3) 

i=0 ^ '' ' 

where fa{ET) = 5{Et) and fi{ET) is the z-th con- 
volution of /i(i?T) 

.U{x) = / dy h{y) h^^{x - y) . (4) 
Jo 

Substituting Eq.|21into Eq. ^ 



and reversing the indices gives a form that is con- the present discussion: 
siderably easier to compute and which is relevant to 

\dETj^CM ^ 
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<7BA 



where 



(7) 



and 



N„- 



"'o(Po) = Y PQWn= ^ Po W« (8) 



is the probabiUty for an inelastic B+A reaction to 
produce no signal on the detector, where wq = hy 
definition. j3| Thus, the detected cross section for a 
B+A reaction is a'j^j^ = aBA{^ — Wo(Po))- 

It is important to emphasize that the acceptance 
of the experimental measurement can be accounted 
for correctly by using the measured N-N cross sec- 
tion in the detector, cr^^, to calculate the probabil- 
ity, Po, that an N-N inelastic collision will produce 
zero signal on the detector: 



Po = 1 



' det 



(9) 



and then taking po into account in the overall detec- 
tor response. Q The properly normalized equation 
for 1 N-N collision on the detector is then given by 
Eq.|2]and the signal for n independent N-N collisions 
on the detector is given by the binomial distribution 
Eq.|3| Thus, the true mean for n independent N — N 
collisions on the detector is: 

{ET)lr' ^ j Et PniEr) dEr = n{ET){l-po) 

which is n times the true mean for 1 N-N collision: 

{EtY"'' = J Et Pi{Et) dEr - (St)(1 - Po) 

(11) 

where (Et) is the mean of the measured Ex dis- 
tribution for 1 detected N-N collision, the reference 
distribution. It is important to contrast Eg. 1101 with 
the mean of the rt-th convolution of the observed 
reference distribution, Eq. 01 



J ETfn{ET)dET^n{Ej 



(12) 



which is n times the observed (Et), as it should 
be, and which differs from the mean of the distribu- 
tion P„{Et) for n independently interacting projec- 
tile nucleons (Eq. I10() by a factor of 1 — po for all 



To summarize, the NCM calculation for a B-l-A 
reaction is given by Eq. ^ which uses the relative 
probabilities, Wn, for n independent N-N collisions 
calculated with the inelastic N-N cross section, <J^^ , 
and takes into account the fraction of inelastic N-N 
collisions, po, which produce no signal on the de- 
tector, by correcting the response function, PniEx), 
Eq. O The result is identical to Eq. El which is the 
sum over convolutions of the uncorrected spectrum 
for i independent detected N-N collisions. The cor- 
rection for the effect of po appears in w'i{po) (Eq.[7J), 
which are thus the relative probabilities for i inde- 
pendent detected N-N collisions. J^J It is important 
call attention to a key point of this section: if instead 
of correcting for the acceptance at the detector level, 
the correction is made by adjusting the input cross 
section of the nuclear geometry calculation , then 
the w'i{po) in Eq.[7|should correspond to a Glauber 
calculation with a^^f 



(1-PoX^ 
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Does using the observed a^J^ in the nuclear 
geometry calculation give the same result? 



In the typical static or Glauber Monte Carlo cal- 
culation 0, Q used to calculate the distribution Wn 
of binary N-N collisions, a collision is defined when 

two nucleons pass within a distance r < \J'o^^Jt^ 

from each other. For any impact parameter h the 
number of collisions for a given nuclear configura- 
tion is calculated and the distribution in the num- 
ber of collisions is obtained by integrating over all 
impact parameters. If the distribution w„ is first 
calculated for a cross section ct;^^, then if a smaller 

"Po)o-i^o?' is relevant. 



cross section, e.g. a^J^ 



1 



the probability for m collisions with cross section 
(T^^, given n collisions with cross section a^^, is 
given by the binomial expansion: 



[n — my. to! 



p^™(l-po)" (13) 



to the extent that the collisions are statistically inde- 
pendent. The new probability distribution w'mipo) 
for m collisions is obtained by summing over all 
n > m, where < m < N^s^x 



N„ 



w'.m{po) = ^ W,iP{m)\r 
n—m 

= E 



[n — to)! to! 



p^™(l-po)'",(l4) 
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or, 

w'rniPo) = (I-PO)" E 7 ^ir^Po"™"^" • 

''^ — ^ in — m ! m! 

n—m ■ ' 

.(15) 

Since the relative probability Wn of n collisions 
with cross section of^^ is normalized: 

A'max 

^ = 1 , (16) 

I 



m— m— 

E ^''-(^'o) - E E 7 — V^po""( 

■'^ n — m ! m! 

m=U 71—1 ^ ^ 



Clearly the w'm(po) in Ea.llSland Eq.[7|are identical. 
This proves that where the probability for a colli- 
sion is proportional to the cross section in question, 
and when the probabilities for individual inelastic 
collisions are independent, that the same result for 
the extreme independent calculation of the Et dis- 
tribution in the NCM model for a B-l-A reaction is 
obtained whether the detector response per collision 
is corrected for the probability po of recording zero 
signal for an inelastic collision, or when the actual 
measured cross section in the detector is used in the 
calculation of the nuclear geometry. This condition 
may also apply for the AQM, and a specific example 
of the AQM with the detected cross section used in 
the nuclear geometry calculation has been given in 
the literature for a — a collisions at the CERN- 
ISR. [l| 

However, the conditions of statistical indepen- 
dence and linearity of the number of elementary 
"collisions" with cross section do not apply for 
the cases of the Wounded Nucleon Model or the 
Wounded Projectile Nucleon model. This is shown 
for the WNM in Fig. ^ where w„ calculated with 
either 42mb or 30mb does not change except for 
the most central collisions, where there is a slight 
difference — this is rather different from the combina- 
torial suppression factor that appears in the correct 
computation of w' nij>n) with 1 — po = 30/42 = 0.714 
(Eqs.EE). Similar logic applies to the WPNM cal- 
culation, where the number of projectile nucleons 
struck by a target nucleus A, with a slightly differ- 
ent p-A cross section, would hardly change for cen- 



it is easy to see that w' m are also normalized, which 
follows by reversing the indices m and n in Eq. 1171 : 



(17) 



-^0)"= E + = l • (18) 

n=l 

I 




n/2=Np„,/2 



FIG. 1: ™„/2 =Prob(iVp<,rt/2) for Au-^Au with cr^ = 
30mb (E802 3]) or 42 mb (PHENIX [ij) The dotted 
hne represents w'„/2 calculated with with ct;^^ = 42mb 
and 1 - po = 30/42 = 0.714. 

tral collisions, and then only for nucleons near the 
periphery. 

One may question even for the number of colli- 
sion model whether the condition for independence 
of collisions applies. In the Glauber model, an in- 
dividual nucleon in a projectile nucleus is defined 
to strike nucleons in the target nucleus when the 
nucleons in the target are found within a distance 

^ < \l o'hiii' / from the line of the projectile tra- 
jectory. In a Glauber Monte Carlo calculation, the 
randomness of assembling nucleons into a target nu- 
cleus with a density 0.16 fm~'^ (following a Woods- 
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n—m. ^ ' 
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Saxon distribution assures the randomness of 

the number of target nucleons in a cyhnder of ra- 
dius 1.16fm (0.98fm) corresponding to cr^^^ — 42mb 
(30mb). Thus the probabihty to find a nucleon in its 
spot of phase space in the cyhnder (i.e. the proba- 
bihty for a binary coUision) for a p-l-A interaction at 
a given impact parameter is random so the coUisions 
are independent. However, for a B-l-A reaction, as 
the impact parameter varies, the total number of 
nucleons involved varies in a correlated manner: the 
more central the interaction, the larger the number 
of nucleons that can make binary collisions. Thus, 
the number of binary collisions is not fully random 
because each of the nucleons participating tends to 
make more or fewer binary collisions depending on 
the impact parameter. Thus, the distribution in 
Wn is not binomial for B+A interactions. Never- 
theless, statistical independence is not unreasonable 
for the change in probability corresponding to the 
change in radius of O.lSfm, which is much smaller 
than the random variation of the number of target 



nucleons along a given projectile trajectory, so the 
binomial distribution for the change in the number 
of coUisions with the change in cross section, Eq.El 
is probably correct. 



III. CONCLUSION 

We have demonstrated that the use of the experi- 
mentally detected cross section in the nuclear geome- 
try (Glauber) calculation is valid in certain classes of 
models but does not work in general. The procedure 
which always works 3] is to perform the nuclear ge- 
ometry calculation with the inelastic N-N cross sec- 
tion and correct the detector response for the prob- 
ability of measuring zero for an inelastic collision or 
other fundamental element of particle production. 
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